Can a supernova be located by its neutrinos? 
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A future core-collapse supernova in our Galaxy will be detected by several neutrino detectors around 
the world. The neutrinos escape from the supernova core over several seconds from the time of 
collapse, unlike the electromagnetic radiation, emitted from the envelope, which is delayed by a 
time of order hours. In addition, the electromagnetic radiation can be obscured by dust in the 
intervening interstellar space. The question therefore arises whether a supernova can be located by 
its neutrinos alone. The early warning of a supernova and its location might allow greatly improved 
astronomical observations. The theme of the present work is a careful and realistic assessment of 
this question, taking into account the statistical significance of the various neutrino signals. Not 
surprisingly, neutrino-electron forward scattering leads to a good determination of the supernova 
direction, even in the presence of the large and nearly isotropic background from other reactions. 
Even with the most pessimistic background assumptions, SuperKamiokande (SK) and the Sudbury 
Neutrino Observatory (SNO) can restrict the supernova direction to be within circles of radius 5° and 
20°, respectively. Other reactions with more events but weaker angular dependence are much less 
useful for locating the supernova. Finally, there is the oft-discussed possibility of triangulation, i.e., 
determination of the supernova direction based on an arrival time delay between different detectors. 
Given the expected statistics we show that, contrary to previous estimates, this technique does not 
allow a good determination of the supernova direction. 
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I. INTRODUCTION 

There has been great interest recently in the question 
of whether or not a supernova can be located by its neu- 
trinos. If so, this may offer an opportunity to give an 
early warning to the astronomical community, so that 
the supernova light curves can be observed from the ear- 
liest possible time. An international supernova early alert 
network has been formed for this purpose, and the details 
of its implementation j^] were the subject of a recent 
workshop ||. (The creation of such a network was also 
discussed in Ref. [f|). One of the primary motivations for 
such a network is to greatly reduce the false signal rate by 
demanding a coincidence between several different detec- 
tors, as detailed in Refs. [0J^]. The second motivation, to 
locate the supernova by the neutrino signal, is the topic 
of this paper. 

An early-warning network is important because super- 
novae are rare, with the estimated core-collapse super- 
nova rate in the Galaxy about 3 times per century [^J. 
The present neutrino detectors can easily observe a su- 
pernova anywhere in the Galaxy or its immediate com- 
panions (e.g., the Magellanic Clouds). Unfortunately, 
the present detectors do not have large enough volumes 
to observe a supernova in even the nearest galaxy (An- 
dromeda, about 700 kpc away). 

Since decisions about how to implement this network 
are being made now, it is of current and necessary prac- 
tical interest to make detailed calculations of what can 
realistically be done. In this paper we carefully examine 
the available techniques for locating the supernova by its 



neutrino signal. The problem has been discussed in gen- 
eral before (see Ref. B, for example). (See also the early 
discussions in Refs. |7 -[)]]). Previous estimates, in par - 
ticular regarding the triangulation problem [p]J^J^, [Loyi]| , 
were rather optimistic. In this paper, we make explicit 
some of the underlying assumptions in these calculations, 
and explain what the fundamental limitations on the pre- 
cision are. Generally speaking, we find that this problem 
is more difficult than had been anticipated. The results 
below were first presented at the above workshop 

There are two types of techniques to locate a supernova 
by its neutrinos. The first class of techniques is based 
on angular distributions of the neutrino reaction prod- 
ucts, which can be correlated with the neutrino direction. 
In this case, a single experiment can independently an- 
nounce a direction and its error. The second method of 
supernova location is based on triangulation using two or 
more widely-separated detectors. This technique would 
require significant and immediate data sharing among 
the different experiments. Our conclusion, that triangu- 
lation is at best very crude, has an important impact on 
the ongoing design decisions for the type of alert network 
and the degree of data sharing. 

Let us note in passing that the supernova might be not 
only located on the sky, i.e., in two dimensions, but that 
its distance can be also reasonably estimated. The num- 
ber of neutrino events N is proportional to the binding 
energy release Eb of the supernova, and of course falls 
off as the distance D squared: 
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where Nq is the number of events at the canonical values 
of Eb and D. The binding energy is thought to be Eb = 
(3. Oil. 5) x 10 53 ergs Q, i.e., a relative precision of 50%. 
The total numbers of events from all reactions are Nq ~ 
10 4 for SK and N ~ 10 3 for SNO (the expected signals in 
SK and SNO are discussed in Refs. Qjjj). Note that iV 
depends on the neutrino spectrum temperatures (which 
can be measured) . The relative errors on the measured iV 
and the calculated No are much smaller than the relative 
error on Eb- Then 



SD 
~D 



ISEb 
2 E R 



(2) 



and thus D can be determined with a relative precision 
of order 25% by any detector with reasonable statistics. 

Finally, it is likely that the next supernova will lie in 
the Galactic plane (including the bulge). However, in 
our opinion, the point of the neutrino measurement is to 
make an unbiased estimate of the supernova location, so 
we do not use this as a constraint. For concreteness, we 
assume that D = 10 kpc (approximately the distance to 
the Galactic center), but in an arbitrary direction. 



II. REACTIONS WITH ANGULAR 
DEPENDENCE 

Neutrinos are detected by their interaction with the 
target material, i.e., its electrons or nuclei. For some 
reactions, the angular distribution of the reaction prod- 
ucts is correlated with the incoming neutrino direction. 
In this section we describe how this angular dependence 
can be used for determination of the supernova direction. 



A. Neutrino-electron scattering: forward peaking 



Neutrino-electron scattering, 



(3) 



occurs for all flavors of neutrinos and antineutrinos, and 
is detected by observing the recoil electrons with kinetic 
energy T above the experimental threshold T m i n . The 
scattering angle is dictated by the kinematics and is given 
by 
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Both SK and SNO hope to have a threshold of order 
T m in = 5 MeV, and so cos a > 0.91. However, after 
integrating over the electron kinetic energy distribution 
for a fixed neutrino energy, and also the neutrino energy 



spectrum, the average value will be larger. We take the 
latter to be of the Fermi-Dirac type with temperature 
T = 3.5, 5 and 8 MeV for v e , D e and v x = v^, z> M , i/ T , v T , 
respectively. Then we obtain (cos a) = 0.98, 0.97, and 
0.98, respectively, with the combined average (cos a) = 
0.98, corresponding to about 11°. 

The angular distribution of the produced electrons is 
narrow, and depends on energy and flavor. However, 
multiple scattering of the electron will smear its Cerenkov 
cone. This washes out the dependence on energy and 
flavor, and one can reasonably model the electrons as 
having a Gaussian smearing from the forward direction, 
with a one-sigma width of 25°, for all energies and flavors. 
This is consistent with the estimate for SNO jl(| and the 
measurement by a LIN AC for SK JlTj . The SK measure- 
ment shows that the angular resolution does depend on 
the electron energy, but the variation is not large. 

Naively, if the one-sigma angular width of this cone is 
da ~ 25°, then the precision with which its center (i.e., 
the average) can be defined given N$ events is 



(5) 



where 8 measures the angle from the best-fit direction 
(i.e., the average). For SK @, N s ~ 320, so the cone 
center could be defined to within about 1.5°. For SNO 
(using both the light and heavy water) ]15|| , Ns — 25, 
so the cone center could be defined to about 5°. These 
results (at least for SK) are widely known (see, for ex- 
ample, Ref. 0). The equivalent error on the cosine is 
<5(cos6<) ~ (50) 2 /2, i.e., 3 x 10~ 4 and 4 x 10~ 3 , respec- 
tively. 

These results neglect the fact that the centroiding is 
to be done in two dimensions, and that this peak sits 
on a large background. It has been claimed that the er- 
ror on centroiding in two dimensions is larger than the 
corresponding error in one dimension by a factor of >/2. 
For r = \J x 2 + y 2 , and uncorrelated errors of equal mag- 
nitude, Sx = 5y = a, simple error propagation gives 
5r = a. Only for correlated errors, e.g., a positive error 
Sx always accompanied by a positive Sy of the same mag- 
nitude, does the factor y2 appear. Centroiding in two 
dimensions is no harder than centroiding in one dimen- 
sion since there are twice as many measurements, i.e., x 
and y for each point. 

However, the nearly isotropic background from all 
other reactions, neglected in previous estimates, is more 
of a concern. Finding the supernova direction becomes 
a question of finding the centroid of a Gaussian peak of 
known width on a known flat background. The centroid- 
ing precision can be determined by a test due to Du- 
Mond |l8| . This result follows from the assumption that 
a known template function with unknown centroid is ad- 
justed until it gives the best least-squares fit to the data, 
i.e., exactly what one would do in practice. This tech- 
nique does not require any arbitrary cuts on which data 
are included in determining the centroid (which would 
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introduce bias). The appropriate two-dimensional gener- 
alization is 
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where L(x, y) = d 2 N / dxdy , the density of events. Note 
that since the error is the same in any direction, and the 
x and y directions are arbitrary, we have considered the 
case in which the centroiding error Sx is only in the x 
direction. One can show that this Sx is exactly the shift 
in the centroid which would increase the total % 2 of the 
fit by 1. (One can also use the Rao-Cramer theorem, 
introduced in a later section, to show that this error Sx 
is the minimum which can be achieved by any technique). 

We use a two-dimensional Gaussian peak (with a total 
of Ns signal events) on a flat background: 
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where R is the ratio of the heights of the flat background 
and the signal (at peak). For Nb background events on 
the whole sphere, 
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where a = 8a ~ 25°. Once R has been specified, we 
treat the problem in a plane instead of on the sphere, 
since a is small. When R = 0, the integral can be done 
analytically, with the expected result of Sx = a / V Ns ■ 
For R > 0, it must be done numerically. We define a 
correction factor C(R) to the naive error by 
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where the full error Sx is determined numerically from 
Eq. (§). The function C(R) is shown in Fig. 1. Note 
that C(R) depends on R alone, so that the same figure 
can be used for both SK and SNO. Empirically, C(R) is 
reasonably fit by C(R) ~ vl~+ AR. This form is moti- 
vated by two constraints: that C(0) = 1, and that for 
R > 1, one can show from Eq. (@) that C(R) ~ V4R. 

For SK and SNO, the number of signal events Ns is 320 
and 25, respectively. As a worst case, we assume that all 
events from other reactions are background events. Then 
for SK, R ~ 3.0 and C(R) ~ 3.7, and for SNO, R ~ 3.8 
and C(R) ~ 4.1. These correction factors may seem 
surprisingly large, but one should note that a two-sigma 
region contains about 2000 and 200 background events 
for SK and SNO, respectively. With the most pessimistic 
background assumptions, the centroiding errors for SK 
and SNO are then about 5° and 20°, respectively. 

It should be possible to reduce this isotropic back- 
ground. In neutrino-electron scattering, the outgoing 
electrons tend to have energies well below the neutrino 
energy. In contrast, in the reaction v e + p — > e + + n, 




FIG. 1. The correction factor to the naive pointing error 
from v + e~ — > v + e~ scattering due to the isotropic back- 
ground is shown. The solid line with points is the numerical 
result. The dashed line is + 4R, an approximation that is 
discussed in the text. 



the outgoing positron carries almost all of the neutrino 
energy. Approximately 2/3 of these background events 
are above 20 MeV, and can be cut with little loss in 
signal. The background in the SNO heavy water will 
depend on the neutron detection technique. Crude esti- 
mates indicate that SK and SNO may each be able to 
attain C(R) ~ 2 - 3. 



B. Neutrino-nucleus reactions: weak angular 
dependence 

In this subsection we consider charged-current reac- 
tions on nuclear targets, i.e., reactions in which only v e 
and v e participate. The reaction with the most events is 
v e + p — > e + + 7i, with ~ 10 4 events expected in SK, and 
~ 400 events expected in the light water of SNO. The 
other relevant reactions are those on deuterons in SNO, 
with ~ 80 events each expected for v e + d — > e + + n + n 
and v e + d — > e~ +p + p fig] . We neglect charged-current 
reactions on the isotopes of oxygen , which also have 
weak asymmetries, but are difficult to separate from more 
dominant reactions. The reactions considered have lep- 
ton angular distributions approximately of the form 
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with, in general, a = a(E u ). We will discuss the magni- 
tude and variation of the coefficient a shortly. We first 
consider how well one could localize the supernova as- 
suming that a is known and constant. In the following, 
we neglect the experimental angular resolution for elec- 
tron and positron directions, as it will be negligible in 
comparison to the pointing errors discussed below. 
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Given a sample of events, one can attempt to find the 
axis defined by the neutrino direction. Along this axis, 
the distribution should be flat in the azimuthal angle cj) 
and should have the form of Eq. ( fio| ) in cos a. Along 
any other axis, the distribution will be a complicated 
function of both the altitude and azimuthal angles. We 
assume that the axis has been found numerically, and 
ask how well the statistics allow the axis to be defined. 
A convenient way to assess that is to define the forward- 
backward asymmetry as 



.4 



FB 



N F -N B 
N F + N B 



(11) 



where Np and Ng are the numbers of events in the for- 
ward and backward hemispheres. The total number of 
events is N = Np + Nb- Note that Apb will assume an 
extremal value A e p^ along the correct neutrino direction. 

The error on Ap B due to the statistical errors on Nf 
and Nb is 
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Using Eq. (J10|) one finds simply, 
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where the error is nearly independent of a for small \a\, 
which is the case for the reactions under consideration. 

In the above, the coordinate system axis was consid- 
ered to be correctly aligned with the neutrino direction. 
Now consider what would happen if the coordinate sys- 
tem were misaligned. While in general, all three Euler 
angles would be needed to specify an arbitrary change in 
the coordinate system, symmetry considerations dictate 
that the computed value of Aps depends only upon one 
- the angle 6 between the true and the supposed neu- 
trino axis. Thus Aps is some function of if the axis is 
misaligned. Using a Legendre expansion, one can show 
that 



Afb iff) = ~ cos( 
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The same expression is obtained by expressing the angu- 
lar distribution in spherical harmonics and considering its 
behavior when acted on by the rotation operator. The 
error on the alignment is then 
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As noted, one would find the best estimate of the neu- 
trino axis numerically and define that direction to be 
cos 6 = 1. Along that axis, the measured asymmetry will 
be A e pg. The above considerations describe the situa- 
tion when the fluctuations do not dominate the value of 
A e pg r ~ a/2 ± 1/y/N. Only in that case one can hope 
to use the angular distribution for pointing, and at the 
same time avoid apparent formal difficulties like an in- 
finite error in the cosine when a — > or a possibility of 
|cos0| > 1 in Eq. @. 

Treating the nuclcons as infinitely heavy, the coefficient 
a in Eq. ( [10| ) is related to the competition of the Fermi 
(no spin flip) and Gamow- Teller (spin flip) parts of the 
matrix element squared: 
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For the v e + p — > e + + n reaction, \Mgt/Mp\ — 1.26 
and thus a ~ —0.1. However, as we have shown else- 
where |Q , due to recoil and weak magnetism corrections 
of order 1/M p , where M p is the proton mass, the coef- 
ficient a varies quite rapidly with neutrino energy and 
changes sign near E v = 15 MeV, becoming positive at 
higher energies (see also Ref. pH). In fact, after aver- 
aging over the v e spectrum, which is taken as before as 
Fermi-Dirac with a temperature T — 5 MeV, we obtain 
(a) ~ +0.08. 

As stated above, for v e + p — * e + + n one expects 
N ~ 10 4 events in SK. This would imply <5(cos0) ~ 0.2 
if a = —0.1 and a similar error if a = +0.08 and the 
temperature is known to be T = 5 MeV (otherwise an 
uncertainty in the temperature would obviously cause an 
additional uncertainty in cos 6). The temperature can 
be measured from the shape of the positron spectrum, 
with precision of order 1% ]l4] , ^5[ . Thus, even though 
the asymmetry parameter a is quite small, the number 
of events is large enough that this technique in SK could 
give a reasonable pointing error. There are also events of 
this type in the light water of SNO. However, N ~ 400, 
so this would imply <5(cos0) ~ 1.0, which is too large to 
be useful. 

The reactions on deuterons in the (M p — > oo) approx- 
imation are pure Gamow- Teller and thus a = — 1/3. We 
will, for the sake of an estimate, assume that for the reac- 
tions on deuterons the energy dependence of a can be ne- 
glected (though see Ref. ^2|). We assume optimistically 
that the signal in SNO of the reactions D e + d — > e + +n+n 
and v e + d — > e~ + p + p, can be combined, so N ~ 160. 
With a = —1/3 this gives 8(cos6) ~ 0.5, which is again 
rather large. 

In the reaction D e +p — > e + + n, the kinematics dictate 
that the the outgoing neutrons have a forward angular 
distribution. If the positions of the positrons and the neu- 
trons can be separately determined, the vector between 
these points can give some information on the neutrino 
direction, at least on a statistical basis ^jj. In fact, this 
effect was observed in the Goesgen p4[. and Chooz |23] 
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reactor experiments. It is not currently possible to detect 
neutrons in SK or the light water in SNO, and we do not 
consider this further. However, this technique may allow 
a scintillator detector to have some pointing ability [E0| . 

The techniques of this subsection will not allow the 
supernova to be located anywhere near as precisely as 
by neutrino-electron scattering. Nevertheless, they may 
provide an independent confirmation of the neutrino- 
electron scattering results, which would increase the con- 
fidence that a real supernova was seen. For example, 
consider the positrons from from v e +p — > e + + n in SK. 
Along the axis determined by neutrino-electron scatter- 
ing, the measured value should be Afb — +0.04, with 
an error of 0.01, a four-sigma effect. This can probably 
be improved somewhat by considering only the highest- 
energy positrons. While the number of events will be 
reduced, the average a will be increased [^o), thus im- 
proving the pointing ability. 

III. TRIANGULATION 

For two detectors separated by a distance d, there will 
be a delay between the arrival times of the neutrino pulse. 
The magnitude of the delay At depends upon the angle 9 
between the supernova direction and the axis connecting 
the two detectors. Given a measured time delay At, the 
unknown angle 9 can be determined: 

cos6»=^. (19) 
d 

The Earth diameter is d » 40 ms. For a typical pair of 
detectors, the time will be somewhat less; for SK and 
SNO, d w 30 ms. The error on the time delay will cause 
an error on the determination of 9: 

5( COS 0) = ^M. (20) 

Thus two detectors define a cone along their axis with 
opening cos 9 and thickness 2 x 5 (cos 9) in which the su- 
pernova can lie. Obviously, in order to have a reason- 
able pointing accuracy from triangulation, one will need 
8 (At) <C d. In what follows we discuss whether an ap- 
propriate time delay can be defined, and what its error 
would likely be. 

Note that this simple error analysis would have to be 
modified near At = ±d, since we must have | cos#| < 1, 
but we ignore this complication. Also, for convenience we 
use cos 9 rather than 9 itself. Naively, 89 = 8(cos 9) / sin 9. 
This is indeed valid for moderate angles, but has spurious 
singula rities at 9 = 0, tt. In fact, for small 89, one has 
89 ~ v /27(cos#) near = 0, it. 

For now we will consider just two detectors, SK and 
SNO, taking events from all reactions in SK and the light 
water in SNO. These are about 10 4 and 400 events, re- 
spectively, mostly v e + p — > e + + n. The effect of mul- 
tiple detectors will be discussed later. Further, we make 



the following assumptions: that the detectors have per- 
fect efficiency at all energies, perfect time resolution and 
synchronization, no dead time, and a negligible time- 
independent background. In practice, these should be 
reasonable assumptions. Therefore, we consider that SK 
and the light water in SNO are ideal detectors, identi- 
cal except for size. The event rates in the two detectors 
should then be identical, except for normalization and 
fluctuations, and a delay At. That is, we are considering 
the best that triangulation could do under any circum- 
stances, limited only by statistical errors. 

The supernova neutrino pulse and hence also the ob- 
servable scattering rate is assumed to have a short rise, 
followed by a relatively slow decline. The total duration 
of the risetime is unknown, but in many supernova mod- 
els it is of order 100 ms p6| . The total duration of the 
decaying phase is much longer, and in most supernova 
models is a few seconds p6[ . However, the observed du- 
ration of SN1987A was clearly longer, about 10 s, and we 
will base our estimates on that. 

If we wish to assume as little as possible about the 
form of the event rate, we could simply consider the shift 
in the average arrival time (t) between the two detectors. 
The expected value is just At. The error on the deter- 
mination of the average arrival time is the width of the 
pulse divided by the square root of the number of events. 
Quite generally, this must be of order 3 s/Vl0 4 ~ 30 
ms for SK and 3/V400 ~ 150 ms for SNO, and there- 
fore not useful for triangulation. One can show that a 
Kolmogorov-Smirnov test for a delay between the SK and 
SNO data leads to the same result for the error. 

In our previous studies of the effect of a v T mass on the 
neutral-current event rate in SK or SNO |l4|,[l5| , we com- 
pared the average arrival time from the neutral-current 
events to the average arrival time from the charged- 
current events. We showed that this allows the detec- 
tion of a delay and the extraction of a v T mass with only 
minimal assumptions about the time dependences of the 
event rates. The results obtained demonstrated the best 
that one could do without assuming specific models for 
the event rates. 

Evidently, the triangulation problem will require more 
detailed assumptions about the time dependence of the 
event rate. As noted above, here we have an additional 
piece of information, that besides the delay At, the event 
rates in SK and the light water in SNO should differ only 
in normalization and fluctuations. As a schematic model, 
we consider an event rate which consists of an exponen- 
tial rise with a short time scale (t% ~ 30 ms), followed 
by slower exponential decay (tq. ~ 3 s). The point of 
transition between the two is labeled to- In Fig- 2, the 
normalized event rate f(t) is shown. We show below how 
the results depend upon the timescales t\ and r 2 . Since 
SK will have many more events than SNO, one could ef- 
fectively measure the event rate in SK and use that to 
replace our assumed form. 
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A. Zero risetime case 



1.0 



We consider first an even simpler model, in which the 
risetime t\ = 0. Then the normalized event rate can be 
written 



/(*) 



1 

Tl 



■ exp 



72 



t>t , 



(21) 



and is zero otherwise. If /(f) is guaranteed to have 
an infinitely sharp edge, as above, one can show (see 
Refs. p7|-p9t) that the best estimator (i.e., the technique 
with the smallest error) of the edge t is 



to 



t T2 



where £i is the measured time of the first event, 
error on the determination of to is 



8(h) 
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Here t 2 /N is simply the spacing between events near the 
peak. For a more general f(t), but still with a sharp edge 
at to, one would simply replace t 2 by l//(to). In fact, 
the shape of /(f) is irrelevant except for its effect on the 
peak rate, i.e., f(t ). So long as f(t) has a sharp edge 
and the right total duration, allowing a more general time 
dependence would therefore not change the results signif- 
icantly. That is, only the long timescale t 2 is important, 
since it determines the event spacing near to- 

The event rate in SNO will consist of N events sam- 
pled from f(t,h), and the event rate in SK will consist of 
N' events sampled from f(t, t' ) (where here we show the 
offset explicitly). The parameters h and t' can be ex- 
tracted from the times of the first events as above. Then 
At = to — t , and its error will be dominated by the er- 
ror on to as extracted at SNO, so that 6 (At) ~ Sto- For 
t 2 = 3 s, this idealized model would allow the offset of 
the edge to be measured to ~ 0.3 ms in SK, and ~ 8 ms 
in SNO. This gives 8 (cos 6) ~ 0.25 at one sigma. 

This technique of using the first event only works if 
there is no time-independent background and if there is 
absolutely no tail of f(t) before h- in either case, the 
fluctuation of a single unwanted event can change the 
extracted delay in an unpredictable way. At the cost of 
an increase in the error, this technique could be made 
robust by looking at the average time of the first few 
events. In all of the other techniques discussed in this 
paper, the role of the time-independent background is 
negligible. 



B. Nonzero risetime case 

The above model has somewhat limited use, since the 
assumption of a zero risetime does not seem to be justi- 
fied. As noted, the supernova models suggest a nonzero 
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FIG. 2. The schematic form of the normalized event rate 
f(t) is shown. To the left of to there is an exponential rise with 
time constant ti. To the right of to there is an exponential 
decay with time constant T2. The tick marks on the t-axis 
are in units of the respective time constants n and r%. For 
clarity of display, we used ri/T2 — 10 -1 in the figure, instead 
of the tl/ts ~ 10~ 2 assumed in the analysis. 



risetime, of order t\ ~ 30 ms, related to the shock prop- 
agation time across the supernova core. If the risetime 
is nonzero, the results of the previous subsection cannot 
be used. The error on to is only given by the spacing 
between events at the peak if the edge is sharp. In this 
subsection, we allow a nonzero risetime and show that 
the results have a qualitatively different dependence on 
the parameters than in the previous case. In addition, 
softening the leading edge will obviously make the trian- 
gulation error larger. 

For the normalized event rate, we take 
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Then f(t) is a normalized probability density function 
built out of two exponentials, and joined continuously at 
t = to- In what follows, we assume that this form of f(t) 
is known to be correct and that t\ and t 2 are known. 

As above, the event rate in SNO will consist of N 
events sampled from /(f, t ), and the event rate in SK 
will consist of N' events sampled from f(t,t' ). Then 
At = t — 1' , and again 8 (At) ~ St , since the SNO error 
dominates. We consider only the statistical error deter- 
mined by the number of counts. Any uncertainties in the 
form of /(f) or its parameters will only increase the error. 
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As noted, we want to determine the minimal error on the 
triangulation. 

This model, while simple, contains the essential 
timescales and an adjustable offset. More general mod- 
els for the event rate must reproduce these timescales in 
order to be physically plausible, and so the final value 
for the error would be close to what is obtained here. To 
define the time delay between the two pulses, we have 
used the offset at which the rising and falling exponen- 
tials are joined. However, in the final result for the error 
on the delay, the particular way in which the offset time 
is defined drops out and the result is therefore quite gen- 
eral. That is, this simple model for the event rate is 
general enough for calculating the statistical error on the 
measured delay. 

These considerations lead to a well-posed statistical 
problem: If TV events are sampled from a known distri- 
bution f(t, to), how well can to be determined? The Rao- 
Cramer theorem |27],|30| provides an answer to this ques- 
tion. This theorem allows one to calculate the minimum 
possible variance on the determination of a parameter 
(here to), by any technique whatsoever. This minimum 
variance can be achieved when all of the data is used 
as "efficiently" as possible, which is frequently possible 
in practice. One requirement of the theorem is that the 
domain of positive probability must be independent of 
the parameter to be determined. This condition is obvi- 
ously not met for a zero risetime, since then the domain is 
(to, oo). For a nonzero risetime, the domain is technically 
(—00,00), independent of to, and so the theorem applies. 
The minimum possible variance on the determination of 
t is: 
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This is the general form for an arbitrary parameter to. 
When to is a translation parameter, i.e., /(t,to) depends 
only on t — to, this reduces to 
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The latter form (the DuMond form) was independently 
derived by another method in Ref. |18| . For the particu- 
lar choice of /(t,to) above, this reduces to 



(St ) 



N x (ai/rl +a 2 /r%) 



For ti <Ct 2 , the minimum error is then 



N 
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(31) 



Note that iVi ~ N(ti/t2) is the the number of events 
in the rising part of the pulse. Since the rise is the 



sharpest feature in f(t), it is unsurprising that it con- 
tains almost all of the information about to. The total 
number of events N is fixed by the supernova binding 
energy release. A change in the total duration of the 
pulse, i.e., ti, would therefore affect the peak event rate 
and hence the fraction of events in the leading edge, i.e., 
Ni/N — ai ~ T1/V2. That is, for N fixed, we are con- 
sidering how a change in the assumed value of ti would 
affect the timing sensitivity; note that T2 appears only 
via the fraction of events in the leading part of the pulse. 
For a more general jit), one would replace ti/t% by this 
fraction computed directly. 

ForSNO, Ni ~ 10- 2 x400 ~ 4, so 5{t Q ) ~ 30 ms/V4 ~ 
15 ms. Since SK has about 25 times more events, the cor- 
responding error would be about 3 ms. Therefore, the 
error on the delay is 8 (At) ~ 15 ms and 8(cosd) ~ 0.50 
at one sigma. We have not yet specified the method for 
extracting to and hence At from the data. That is ex- 
actly the point of the Rao-Cramer theorem - that one 
can determine the minimum possible error without hav- 
ing to try all possible methods. A possible technique 
which should come close to achieving this minimal error 
is discussed below. 

For the two cases, zero and nonzero risetime, we used 
different mathematical techniques which were applicable 
only in one case or the other. This may seem like an artifi- 
cial distinction, and that these two cases do not naturally 
limit to each other. In particular, it may seem incompat- 
ible that in the first case the error ~ l/N, while in the 
second the error ~ 1/yiV. Further, one obviously cannot 
take t 1 — ► in the results of this subsection. 

However, there is physically a natural joining of the two 
results. For the case of t\ — 0, we found that 8 (to) ~ 8 
ms for SNO. For the case of t\ > 0, we chose t\ = 30 
ms and found Ni ~ 4 and 8 (to) — 15 ms for SNO. If we 
reduce r\ by a factor 4, so that t\ — 7.5 ms, then this 
error is reduced by a factor 2 so that 8 (to) ~ 7.5 ms. But 
now there is only Ni ~ 1 event in the rising part of the 
pulse. At this point, the difference between t\ = and 
Ti < 7.5 ms becomes difficult to distinguish, i.e., the edge 
appears sharp. That is, the two techniques give the same 
numerical result for the error at the boundary between 
the two cases. 



C. What will the event rate really look like? 

We considered two simple models for what the event 
rates might look like. Those models were of course crude, 
yet they illustrate how the different timescales affect the 
final results. The short rise (~ Ti) at the beginning of the 
pulse is a prominent feature, and it provides most of the 
timing information. The long decay (~ T2) is important 
principally through how it affects the normalization, i.e., 
the number of events before or near the maximum. 

The risetime, set by the timescale tj., may be smaller 
than suggested in Refs. EG]. For example, in some of 
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the models considered in Refs. |3lf] , the duration of the 
rise does appear to be shorter than considered here. In 
some cases, the v e luminosity rises nearly instantaneously 
to a given value, and then more slowly to a peak value 
which is several times larger. In this case, a zero-risetime 
analysis may be appropriate, and the error on to depends 
on the spacing between events. However, one must not 
use the peak event rate, but rather the event rate at the 
point where the rise is no longer instantaneous. This will 
give an error several times larger than the zero-risetime 
case used above, where the instantaneous portion rose all 
the way to the peak. 

The neutrino pulse duration, set by the timescale T2, 
may also be smaller than assumed here, and our results 
can easily be scaled appropriately. Note that in the zero- 
risetime case, the error ~ ti , and in the nonzero- risetime 
case, the error is ~ \fr-i- Thus in the latter case even 
if Ta = 1 s, the triangulation error would improve by 
only a factor y/3, and so 6(cos 9) ~ 0.3, still rather large. 
However, one has to keep in mind that our choice T2 = 3 
s was motivated by the SN1987A observation that 25% 
of the events arrived at least 5 seconds after the start of 
the pulse. While the statistics were poor, a timescale of 
T2 = 1 s or smaller seems to be unlikely. It may also 
be that the decay of the neutrino pulse is characterized 
by two timescales - a quick drop, followed by a slower 
decline, e.g., a simplification of the gradually decreasing 
timescale considered in Rcf. |j32| . Even so, if the duration 
of SN1987A is reproduced, there would be little change 
in the error on the delay. 

The event rate may also have a much more complex 
structure than assumed here. For example, there could 
be oscillations or other sharp features which could be 
used to define a time from which to measure the delay. 
However, since the most prominent possible sharp fea- 
ture, a zero-risetime edge, is not enough for a successful 
triangulation, any such features should be less significant. 
Moreover, one should not forget that the whole point of 
the supernova early alert network is to use an automated 
analysis to determine the direction from the data. For the 
result to be available essentially immediately, the analy- 
sis should assume as little as possible about the shape of 
the pulse. That was part of our motivation to consider a 
simple model characterized only by timescales which are 
reasonably well-known. 

Note also that we neglected the events in the prompt 
burst of electron neutrinos. The expected number of 
events is very small (for example, Ref. J3j| has N ~ 1 
in SNO), although there is considerable variation in the 
predictions for the number of events and the duration of 
this burst. It seems unlikely that the timing error would 
be small enough to be useful, or that the predictions are 
robust enough to allow an automated analysis. 

In any case, the scenario for attempted triangulation 
could be the following: the detector with the largest 
statistics, e.g., SK, could be used to determine the frac- 
tion Ni/N of events in the leading edge and the duration 
of the risetime t\ (both of these will have some errors, un- 



like what we assumed in using the Rao-Cramer theorem). 
This fraction would then specify which of the SNO events 
were to be considered as coming from the leading edge. 
The delay could then be determined by the time differ- 
ence of the average arrival times of leading edge events 
in SK and SNO. The error would then be ~ n/y/N^, 
where Ni is the number of events in SNO coming from 
the rising edge. Thus this technique might approximately 
attain the Rao-Cramer lower bound on the error. Al- 
ternatively if SK determines the zero-risetime model is 
applicable, the delay could be extracted from the time 
difference of first events in each detector. In this case, 
the error is again determined by SNO, and would be the 
event spacing near the peak. 

D. Comparison to other work 

In Ref. jnj , the triangulation error for a zero-risetime 
pulse using SK and SNO was also considered, with a 
final result of 1.3 ms, to be compared with our result of 8 
ms. The difference is due to different input parameters. 
We assumed 400 events in the light water of SNO, an 
exponential decay of the event rate with T2 = 3 s, and a 
distance of D = 10 kpc. In Ref. |hJ, it was assumed that 
all flavors and all reactions could be combined (we argue 
against this below), for 10 3 events in total (at a distance 
of 10 kpc), with half in the first 1 s. The greater number 
of events and the shorter assumed duration of the pulse 
make the event rate at peak 500 s~ 1 instead of 133 s _1 . 
In addition, the final error was scaled to a distance of 8 
kpc. After correcting for these differences, the results are 
in agreement. The same holds for the results in Ref. ||, 
where again a peak rate of 500 s _1 for SNO and a zero 
risetime were assumed. 

In Ref. |y]], a different technique was proposed, which 
does not explicitly specify whether or not the risetime is 
nonzero. The proposed technique begins by construct- 
ing the cumulative distributions (this function increases 
by a step of 1/N at each event time, and is discrete but 
not binned) for SK and SNO. At least for the light water 
events, these functions should be the same up to fluctu- 
ations and a possible delay. The proposal is to make a 
simple low-order fit to the beginning of each cumulative 
distribution and to extract the delay from the difference 
of the intercepts. A preliminary delay error of order 5 ms 
was presented (for D — 10 kpc). However, the fit to the 
cumulative distribution function does not yet take into 
account the fact that the errors on successive steps are 
highly correlated (because most of the data at a given 
step are the data from the previous step). Taking this 
into account will increase the error. In any case, the error 
from this proposed method cannot be smaller than that 
from the time of the first event (zero risetime case) or the 
Rao-Cramer result (nonzero risetime case), whichever is 
appropriate. 
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IV. CONCLUSIONS AND DISCUSSION 

The final uncertainties, calculated for a canonical su- 
pernova 10 kpc away and with a total energy release of 
3 x 10 53 ergs, are summarized in Table I. These are all 
one sigma errors, though larger confidence regions may 
be necessary for making a search for a supernova. 

Neutrino-electron scattering has the best pointing pre- 
cision. Moreover, the calculated precision is largely inde- 
pendent of assumptions about the supernova model. In 
particular, it is totally independent of the time depen- 
dence of the event rate. The isotropic background from 
other reactions degrades the precision somewhat, but it 
is still the most precise technique. 

The angular asymmetry of positrons from the v e +p — > 
e + +n reaction, even when combined with the high statis- 
tics of SK, does not give a comparably small pointing 
error. It makes it possible, however, to check that the 
signal is indeed coming from the right direction. The 
angular distributions from the charged-current deuteron 
reactions are even weaker. 

Under realistic assumptions about the numbers of 
events and the timescales, triangulation with SK and 
SNO appears to be very difficult, if not impossible. Other 
tests for the time delay can be considered. However, 
for either a zero risetime or a nonzero risetime, we have 
shown in the previous section what the smallest possible 
errors are. For fixed inputs, there is simply no way to do 
better. 

But, so far in this paper, all of the concrete results 
were based on using just two detectors, SK and SNO, and 
taking all events in the light water (these are dominated 
by the charged-current signal v e + p — > e + + n). Can 
the pointing accuracy, in particular for the triangulation 
technique, improve if other existing or planned detectors 
are used? 

First, we stress again that for the method to suc- 
ceed the signals in different detectors must differ only 
in the normalization and fluctuations, and a possible de- 
lay. For example, that precludes including the neutral 
current events (dominated by f M , v T , and their antipar- 
ticles) from the heavy water portion of SNO. That is 
because one cannot guarantee that the time dependence 
of the scattering rate for these events is the same as for 
the events in the light water. In fact, at the crucial early 
times, the supernova models suggest that there are dif- 
ferences among the flavors. Since the time dependences 
of the luminosities and temperatures are not known to 
the needed high precision, these differences cannot be 
corrected for. 

There are several detectors, existing or under con- 
struction, that will observe a few to several hundred 
D e +p — > e + +n events. In particular the existing MACRO 
and LVD detectors in Gran Sasso, and Borexino and 
KamLAND detectors under construction will be clearly 
able to "see" a Galactic supernova (see Refs. Q-g^, re- 
spectively), and can undoubtedly contribute to the false 



signal elimination in the planned early alert network. 
However, since these detectors are based on scintillation 
rather than on Cerenkov light, it is not a priori clear 
that the basic requirement of the similarity of response 
is indeed fulfilled. And even if that difficulty could be 
overcome, the numbers of events in those detectors will 
be comparable to what is expected for SNO. Thus our 
estimate of the uncertainty associated with triangulation 
in the last line of Table I is valid for them also. With 
three detectors, there is in principle an improvement in 
the pointing from triangulation. However, that improve- 
ment is minimal if the pointing error from any two de- 
tectors is of order half of the sky. 

The AMANDA detector (or its successor) can perhaps 
observe a supernova via the v e + p — ► e + + n reaction 
in a very large target volume fl38| . However, the prin- 
ciple used for supernova detection will be a fluctuation 
in the (large) background rate. Over an interval of a 
few seconds, the background events dominate the signal 
events, N B > N s , but N s > \/Nb. The actual su- 
pernova events can only be distinguished in a statistical 
sense. Under these circumstances, it will not be possible 
to map out the event rate well enough to make a precise 
measurement of to or some other appropriate time. 

What would it take to make triangulation viable? As a 
simple example, we consider SK and a hypothetical sec- 
ond detector, also with ~ 10 4 events at 10 kpc, and a 
separation of 30 ms between the two. The second detec- 
tor might be very similar to SK, or it might be primarily 
sensitive to neutral-current reactions (in the latter case, 
we assume, despite the strong cautions above, that the 
charged-current and neutral-current event rates can be 
directly compared for timing purposes). For the event 
rate assumed in our main analysis, each detector would 
have a timing error of 6 (to) ~ 3 ms, so that the triangula- 
tion pointing error would be S(cos9) ~ 0.15. Results for 
the more general case of combining several detectors with 
different timing errors have been given in Refs. PJToj. 

The pointing error from the angular distributions al- 
ways scales with 1/y/N. Under the assumption that 
the event rate risetime is nonzero |2{J, the triangulation 
pointing error also scales with l/yiV. Since N ~ 1/-D 2 , 
all of the errors scale linearly with the distance D. The 
triangulation measurement may then become feasible if 
the distance to the supernova is significantly less than 
10 kpc. However, all of the other techniques improve 
by the same factor. (However, if the risetime were van- 
ishing | pl| , then the triangulation error would scale as 
1/N ^LP; see the discussion above.) 

Thus our analysis shows that a Galactic supernova can 
indeed be located by its neutrino signal, and that among 
the possible methods, the best technique by a large mar- 
gin is neutrino-electron scattering in a water Cerenkov 
detector. Currently, either SK or SNO can separately 
make this measurement. In the above, we considered 
the directional information from the neutrino data alone. 
The operators of the alert network or the astronomers 







themselves can of course combine these results with a 
Galactic model of where a supernova is likely to be. Our 
results indicating that triangulation will be very difficult 
do not mean that the data sharing among different detec- 
tors is not worthwhile. Only a coincidence of two or more 
detectors can eliminate false alarms, and be the basis of 
a reliable early alert system. 
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TABLE I. One-sigma errors on how well the direction to the supernova is defined by various techniques, at D = 10 kpc. The 
other parameters used are noted in the text. For neutrino-electron scattering, the most pessimistic background assumptions 
were used. 



technique 


error 


v + e~ forward scattering (SK) 


56 ~ 5°, S(cos8) ~ 4 x lO' 6 


v + e~ forward scattering (SNO) 


56 ~ 20°, <5(cos0) ~ 6 x 10" 2 


v e + p angular distribution (SK) 


5(cos6) ~ 0.2 


v e + p angular distribution (SNO) 


5(cosd) ~ 1.0 


v e + d,i> e + d angular distributions (SNO) 


5(cosd) ~ 0.5 


triangulation (SK and SNO) 


5(cos6) ~ 0.5 
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